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^ | Abstract 

We use a method, inspired by Pohozeav's work, to study asymptotic behaviors 
of non-variational elliptic systems in dimension n > 3. As an application, we prove 
removal of an apparent singularity in a ball and uniqueness of the entire solution. 
■ All results apply to changing sign solutions. 

A classical work by Gidas and Spruck [6] asserts that any nonnegative solution to 
Am + \u\ a ~ 2 u = in R™ with 2 < a < (subcritical case) is trivial. For a = 
Caffarelli, Gidas and Spruck [I] proved that any nonnegative solution in M. n is of the 
form u = (a + fel^l 2 ) -21 ^ - , where a, b are constants. Such problem for elliptic systems are 
also studied, for example, in the studies of Lane-Emden type systems; see Zou [T7] and 
Polacik, Quittner and Souplet pi)] and the references therein. 

On the other hand, the behaviors of changing sign solutions are more delicate. For 
example, there exists a sequence of changing sign solutions to Au + \u\ a ~ 2 u = in IR n with 
^ < a < ~^2> see El- I n this paper, we study under what circumstances a solution to an 
elliptic system in an exterior domain is asymptotic to |x|~ < -™~ 2 ' ) at the infinity. Such decay 
is optimal in the sense that the infinity is a regular point in the inverted coordinates. It 



is known [7] that there exist solutions to Am + u a ~ x = in M. n which decay slower than 
|x|~( n ~ 2 ). Thus, a suitable integrability condition is necessary to exclude such case. 

While the study of changing sign solutions to elliptic systems is interesting by itself, the 
problem is well-motivated by differential geometry. For example, the decay of curvature 
tensors was studied for Yang- Mills fields [16] , Einstein metrics |T] and other generalizations 
|15j . [5], just to name a few. A typical system for the problem is of the form 

A(Rm) ijk i = Q ijk i(Rm,Rm), 

where Rm is the Riemannian curvature tensor and Q is a quadratic in Rm. A natural 
geometric assumption is that \Rm\ is in L2. Therefore, \Rm\ vanishes at infinity and the 
problem is to find out the decay rate. The study of geometrical systems is more subtle 
as (Rm)ijki satisfies an extra relation, the Bianchi identity, and the underlying spaces are 
not Euclidean. 
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The technique we use in this paper is based on the method developed in [5] on asymp- 
totically flat manifolds, where a special geometric setting is considered. In this paper, we 
study general non-variational elliptic systems of the reaction-diffusion type. Our result 
applies to changing sign solutions and includes the supercritical case (i.e., Au + Cu^ 1 = 
with a > -^j, where C is a constant). 

Let V = (Vi, ■ • • , V m ) and f l : M. m — > R. Consider the system of equations 

TO 

^AjAV^fiV), (1) 
j'=i 

where A is a constant positive definite symmetric matrix and % = 1, ■ ■ ■ , m. The system (JTJ 
describes the steady states of the reaction-diffusion systems. The matrix Ay represents 
the diffusion rate and f l (V) is the reaction term. Note that f l (V) or Vi may have no sign. 
We assume the following structure conditions 

(Al) \f(V)\ <C\V\«, 

(A2) \Vf l (V)\ < C\V\ q -\ 

Let K be a compact subset in R n . 

Theorem 1. Let q > and p = |(g — 1). Suppose that f l satisfies (Al) and (A2). 
Let V G L p (R n \ K) be a™ solution to (QJ) in W n \ K. Then \V\ = O(|o;|^^~ 2 >) and | W| = 
0(|x|^ ( ' n_1 - ) ) at the infinity. 

An immediate consequence is a result on singularity removal for affine invariant equa- 
tions. For scalar equations, the problem was studied in [6], [3], [1]. 
Let B\ be the unit ball centered at the origin. 

Corollary 1. Suppose f l are homogeneous functions of degree Let V G L : ^(B\) be 
a solution to (QP in B\ \ {0}. Then V can be extended to a smooth solution to (QP in Bi. 

By performing a linear transformation Wi = J2j AijVj, the system (CQ) can be reduced 
to an equation of the diagonal form AW = fiW). The assumptions (Al)-(A2) and other 
conditions on V or /* equivalently hold for W and /. Therefore, for Theorem [T] and 
Corollary [TJ we may assume without loss of generality the equation is of the diagonal 
form. 

We turn to study the uniqueness of entire solutions for variational systems. Let P(V) 
be a homogeneous function of degree q + 1. Suppose that /' = in ([rj. Note that 
PI) can not be reduced to a variational system of the form AW = f{W) now; it can 
still be reduced to one with diagonal matrix Ajj. For scalar equations, there is a large 
literature on the uniqueness problem; e.g. Gidas and Spruck [6], Bidaut-Veron [2] and 
Serrin and Zou [2|; see also [12] and the references therein. For systems, when P(V) < 
and q > ^| (supercritical case), the problem was studied by Pucci and Serrin [TT] under 
some asymptotic assumption of V. Their result also holds for non-homogeneous function 
P (and more general P(x, V, VV)) satisfying some inequality. 
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Theorem 2. Let q > q ^ ^| and p = ~(q — 1). Suppose P(V) is a homogeneous 
function of degree q + 1. Let V G L p (U. n ) be a solution to (T7J) in R™ with f l = J^-. Then 
V = 0. 

We give the outline of proofs. To fix notations, we denote by dx the volume element 
in M. n and by dS the area element of a hypersurface in R n . Let B r (x) and S r (x) be the 
ball of radius r and sphere of radius r centered at x, respectively. When x is at the origin, 
we simply denote by B r and S r . 

The idea of proof of Theorem [1] is to compare the size of f Rn \ B \ W\ 2 dx (as a function 
of r) to its derivative — j s \W\ 2 dS. Then by ordinary differential inequality lemma, we 
get the optimal decay of jVV^I and as a consequence the decay of \V\. In order to relate 
above two integrands, we use some version of Pohozaev's identity for non-variational sys- 
tems. Pohozaev's ingenious idea [D] is to use a conformal killing field to prove uniqueness 
in a star-shaped domain. This idea was generalized nicely by Pucci and Serrin [11] to 
general variational systems. Our use of the identity is different from the original one. We 
apply the identity to an unbounded domain (the complement of a large ball) and use only 
the size of |/*|. Therefore, our method can be applied to non-variational systems. 

The proof of Theorem [2] is a combination of Theorem [T] and Pohozaev's original idea. 
Since the solution decays fast enough at infinity, no terms from infinity contribute to the 
main integrand. We use the identity differently such that we obtain the uniqueness also 
in the subcritical case, in contrast to the problem in star-shaped regions where one has 
to restrict to the supercritical case. 

Finally, we show that the assumptions in above Theorems are sharp. 

n+2 

Example 1: Consider the equation Au + u™- 2 = in IR n . By [3], nonnegative solutions 
are of the form u = (a + 6|a;| 2 ) _Ii 2 _ . Therefore, u decays as |a;| _ ^ n_2 ' ) at the infinity. This 
example shows that in Theorem [21 the assumption q ^ ^| is necessary. Consider instead 
the equation in B\ \ {0}. There exists a nonnegative radial singular solution with the blow 
up rate \x\ 2~ near the origin. Therefore, in Corollary [TJ the condition V G L~ =5 (5i) is 
sharp. 

Example 2: Consider Au+u q = in MJ 1 . For q > there exists a solution asymptotic 

to |x| - «- T at the infinity; see [7]. Hence, in Theorem [TJ the conditions q = '^p- and V £ L p 
are sharp. Moreover, in Theorem [2} the condition g = ^is also sharp. 

1 Preliminaries 

We collect some standard results in elliptic regularity theory and ordinary differential 
equations. Lemma [TH follow by an argument similar to pQ, section 4. 

Let C s be the Sobolev constant and 7 = Suppose that the nonnegative function 
u G C ' 1 satisfies Am + Cou q > weakly in the sense that 

j (-(V«, V0) + C u q (f))dx > 
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for all < G C£°. Let ip > be a function with compact support and s > 1. Then by 
Cauchy inequality, 



J if 2 u q+s - l dx > Co 1 J p g g2 ^ \<pVu*\ 2 + ^<fu%(V<f, Vui))dx 
By Sobolev inequality, we have 



1 \Vf\ 2 u s )dx. 
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( J {f 2 uydx)^ <C J i^Yf 2 ^ 3 " 1 + (1 + ^^)|V<p|V)<fc, (2) 

where C = C(n,C s , Co). 

In Lemma [TJ131 u is a C 0,1 function. 

Lemma 1. Let p > 1 and g = 22±™ Suppose that the nonnegative function u G L p (B r ) 
satisfies Au + Co it 9 > weakly in B r . Then there exists e > such that if J B u p dx < e, 

71 V 

then sup Br u < Cr~ p ||w|| j>(B r ), where C = C(n,p, C s , Co). 
Proof. Let s = p in (T5]) . Then 

( f(ip 2 u p ydx)^ < C I \u q -\f 2 u p ) + \Vf\ 2 u p )dx 



< C( u p dx)~( / (if 2 u p ydx)- +C \Vif\ 2 u p dx. 

J {supp ip} J J 

if is chosen to be a cutoff function such that ip = 1 in B r / 2 and if = outside B r with 
\Vf \ < Cr~ l . We get 

( f u pl dx)^ < f u p dx. 

J B r/2 r J B r 

Choose a sequence = (2 _1 + 2~ k )r. Apply (and rescale) the above inequality for B rk 
and B r with p k = pry . By Moser iteration, we have sup Br u < Cr p \\u\\LP(B r )- D 

Lemma 2. Let p > and q = 2£±2. Suppose that the nonnegative function u G 
L p (M. n \ B r ) satisfies Au + Cou q > weakly in M. n \ B r . Then there exists e > such that 
if J" Rn \ B u p < e, then u = 0(|x|~ A ) for all A < n — 2 as \x\ — > oo. 

Proof. By Lemma [H it = 0(|x|" n/p ). Let s = p^ > 1 in Then 

( / ip^vPdx)^ <C{( u p dx)^{ j (ip 2 u pn ^ydx)^ + C I \Vf\ 2 u p ^dx. 

J J {supp ip} J J 
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ip is chosen to be a cutoff function such that tp = 1 in B r > \ B 2t and ip = outside -E?2r' \ B r 
with \Vip\ < C(l/r + 1/r'). Let r' oo. Then 

( f V ^u p dx)^ <C{( |V<^| n cia;)^( / u p dx)^. 

J J J {supp Vip} 



And thus 



vPdx)^ < C{ / u p dx)^. 

l \B 2 r JB 2 r\Br 

This gives f Rn \ B u p = 0{r~ 5 ) for some small 5 > 0. Therefore, by Lemma [IJ u — 

0(\x\~p~p). Let A = sup{A : u = 0(|a;|~ A )}. By iteration and a contradiction argu- 
ment, we get that X = n — 2. □ 

Suppose that h > is a C° function. The nonnegative function u G C 0,1 satisfies 
An + CqHu > weakly if 

y (-(V«, V0) + C hu(p)dx > 

for all < G C °°. 

Lemma 3. Let p > 1 and t > ~. Suppose that the nonnegative function h G L t (B r ) satis- 
fies J B h l dx < 2t- n ■ Suppose also that the nonnegative function u G L p (B r ) satisfies Au+ 

1 n 

C hu > weakly in B r . Then sup Br u < Cr p u where C = C(n,p, C S) C , Ci). 

Proof. The proof is by standard Moser iteration. See Morrey [8]. □ 
The following is a basic result in ordinary differential equations; see [5] . 

Lemma 4. Suppose that f(r) > satisfies f(r) < — ~/'(r) + C 2 r~ h /or some a, 6 > 0. 

(a) a b. Then there exists a constant C3 s?/ca t/iat 

/(r) < C 3 r- Q + ^-r- 6 . 

a — 

Therefore, f(r) = 0(r _mm ^ a '^) as r — > 00. 

(b) a = b. Then there exists a constant C3 such that 

f(r) < C 3 r- a + a C 2 r~ a In r. 
Therefore, f(r) = 0(r~ a In r) as r — > 00. 
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2 Proof of Theorem H 

Proof of TheoremUl We first derive a version of Pohozaev's identity for non- variational 
systems. Let Q be a domain in M. n and N be the unit outer normal on dQ. We will perform 
integration by parts repeatedly. 

/ ^2f k (V)x l D l V k dx= I J2 A kAV^i D iV k dx 
J Q k,i J Q j,k,i 

= f - V] AuDiVjDifaDiVjdx + f J2 AkjDiVjXtDiVkNidS 

= / (- 2 AkjDtVjDiVk - J2 DtAyDiVjDiVktydx + [ Yl A^D^x^V^dS 

ij,k i,j,k,l ~ ^ m i,j,k,l 

= (|-1) / ^A^DiVjDiVkdx- ! l -Y,A kj D i V j D iVk x l N l dS 

i,j,k dQ i,j,k,l 

+ f V A^DiVjXiD^idS. (3) 
Jdn i,j,k,i 

It is worth mentioning that X\Di is a conformal killing field in R n . 

As we explained in the introduction, without loss of generality we may assume the 
equation is of the diagonal form, i.e., 

AVi = f\V). (4) 

Note that \V\ and | W| are C 0,1 functions. By fll) and (Al)- (A2), we have 

A|V| > -C\V\ q ; 
A|W| > -C\V\ q - 1 \VV\ 

weakly. Since V G L p (M. n \K), there exists a large number R such that f Rn \ BR \ V\ p dx < e, 
where e is as in Lemma [TJ Applying Lemma [T] to BJxq) where \xq\ > 2r > 2R, we get 

n 

\V\ = OQx\~p). 

Case 1. If ^±2 < q < ^ (or equivalently, 1 < p < ^), then - > n — 2. By Lemma[TJ 
we have |V| = 0(|x| -n / p ). Let ip be a cutoff function such that </? = 1 in B r and ip = 
outside i?2r with |Vy?| < Cr -1 . Applying y9Vi to (j4j) and integrating gives 

/ \VV\ 2 dx < C [ \V\ q+1 dx + ^ / |^| 2 rfx = 0(r n ~ 2 -^) < 0(r- n+2 ), 

where \xq\ > 2r ^> 1. By Lemma [3] with /i = l^l 9-1 , we obtain |VV| = 0(|x| _ ^ n_1 ^) and 
thus \V\ = 0(|a;|-( n - 2 )). 

Case 2. If < q (or equivalently p > r?o), by Lemma [21 |V| = 0(|x| _A ) for all 
A < n — 2. Therefore, 

/ \VV\ 2 dx<C [ \V\ q+1 dx + ^ [ \V\ 2 dx = 0(r n - 2 - 2X ), 

JB r {x ) JB2r{xQ) T JB 2 r{x ) 
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where \xo\ > 2r ^> 1. Moreover, \V\ G L p ' for all p' > - J ^. Choose p' < p close to 
Hence, q > 2p +n . We can then find q' > § such that 



_B r (x ) 



y2q'—n 



where |x | > 2r ^> l.This is possible because A is close to n — 2. By Lemma [3], we obtain 

sup |w| < ^||vy|U 2 ( Br (, )) = oir-"- 1 ), 

where \xq\ > 2r 3> 1. 

Let ft = £ fl \ 5 r in ©. We have 

/ Y d f\V)x l D l V k dx = &-l) [ \VV\ 2 dx- [ lY^lVVfxMdS 
Jo. kl 1 Jn Jan 1 l 

+ [ yZDiVjxMVjNidS. (5) 



Note that 



lim / R\VV\ 2 dS = lim 0(iT 2A ~ 2+n ) = 0. 



Let R — > oo in (jSJ). Then there is no boundary term coming from the infinity. We can 
choose Q = M. n \ B r . The boundary terms only occur on S r . On dfl, N = —-. Hence, 



l \B r 



V/ fc (y)x,D,V fc dz = (~-l) / |VVfdx+ / ^|W| 2 dS-r / 

fe z ^ JR n \B r J S r * J S 



- v 2 



/ |W| 2 dx + 




\VV\ 


2 dS 


/R n \_B r 








/ |Vr| 2 c/x- 


j r 


|W| 


2 dS. 


/R n \B r 









Let G(r) := J M „^ B |Vy| 2 o?x. Since G'(r) = — f s \W\ 2 dS, the previous formula becomes 

G(r) < ^— G'(r) + / £ fiV^D^dx. 

n — z n — z j K n\ Br fc ; 

The key idea is to compare the size of G(r) to that of G'(r). The coefficient in front 
of G'(r) plays an important role. Here is the only place we use the condition of |/*|. We 
have 

/ y2f k {V)xiDiV k dx< [ \V\ q \x\\\7V\dx = 0(r- x{q+1)+n ). 

J~SL n \B r k I JR n \B r 



Thus, 



G(r) < — — -G\r) + Cr-^ q+l)+n . 
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Since q > an d A is close to n — 2, we have A(g + 1) — n > n — 2. By Lemma HJ this 
implies G(r) = 0(r~( n ~ 2 >). By Sobolev inequality, we get 

2n 

|\/| — rf x = 0(r" n ). 

B2r\Br 

Finally, by Lemma □ and El we obtain |V| = 0(\x\-^ n - 2 ^) and |VV| = 0{\x\~ { - n - v >). □ 



3 Proofs of Corollary HJ and Theorem [2 

Proof of Corollary [H Since the equation is invariant under inversion, we transform the 
solution to W 1 \ Bi and apply Theorem [TJ 

Let y = -Ss. Define Ui(y) = . 1 2 V(-A-). This is called the Kelvin transform with the 
property that 

This can also be viewed as the conformal change formula of the conformal Laplacian with 
zero scalar curvature. Therefore, Ui(y) satisfies 

J^AiAyU.iy) = lf*(\y\»-*U{v)) = f{U{y)) 

in R™ \ Bi, where we use that f l is homogeneous of degree Moreover, 

\U\&dy = [ (\V\\y\- n+2 )^dy= [ (| V| \x\ n - 2 )^ \x\- 2n dx 



2h 



I V| n ~ 2 dx < +00. 



Bi\{0} 



Now we apply Theorem [T] with p = and q = We get \U\ = 0(\y\ ( n 2 ^) and 
\VU\ = 0(\y\-^). Hence, |V| = 0(1) and |VV| = OUx^ 1 ). As a result, V G L°° {B x ) 
and VV G L p (5i) for all p < n. 

We show that V is a weak solution to flTJ) in B\. Let G i/g (£>i, IR m ). Let be a 

compactly supported function in B\ \ {0} such that rjk — > 1 a.e. in £>i and ||?7fc|U«(Bi) — ► 
as — > oo. (Such functions were used by Serrin [13].) Then 

/ rjk^AijDiipjDiVjdx = / -^/*(V)<p»77fcda;- / ^ A?7feAi^i A Vda;. 

The last term can be estimated as follows. 

I / y~] Di7] h A ij Lp j DiV i dx\ < CH^H^j^ ||VV||£a(Bi)lhfc|U»(Bi) < C\\r] k \\ Ln{Bl) -> 



as k — > oo. Hence, in the limit 

Thus, V is a weak solution in B\. It follows by elliptic regularity that V G C°°(Bi). □ 
Proof of Theorem^ Let Q = B R . Therefore, N = § in fl3J. We get 




On the other hand, we also have 

/ ( ff +l)P(V)dx= / ^Ii7^ x 
Vs fl Jb r V ^ 




where we use the Euler formula for homogeneous functions. 
Case 1. n > 4. 

By Theorem [TJ when R — > oo, ([S]) becomes 
/ -nP(V)tte = VA fci A^A^x + 0(i?-( n - 2 )) + 0(i?^ +1 ^- 2 ) +n ) 

= (?-!)/ ^iA^AVfc^ + o(l), 

where we use conditions on p, g and n > 4 to get (g + l)(n — 2) — n > 0. Similarly, ([7]) 
gives 

/ (g + l)P(V)dx = - [ S'AuDiVjDiVkdx + 0(iT ( "- 2) ). 
Jb r Jb r i)j>k 
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Combining these two formulas and noting that q + 1 7^ 5 we h na Hy arrive at 



/ Y^A k] D i V 3 D i V k dx = o{l). 



i,j,k 

Since A is positive definite, we have | VV| = and hence V = 0. 
Case 2. n = 3. 

Note that sup \ V\ < ,2/ p \\V\\lp- Combining this fact with Theorem CD, we have \V\ = 
0(\x\~ x ), where A = max{l, - }. Therefore, 

3 2v 6 

X(q + 1) - 3 > maxig - 2, ~(q + 1) - 3} > max{-l + — , -1 + -} > 0. 

p 3 p 

Then (E]) becomes 

/ -3P(V)dx = d-l) [ J2A kj D l V 3 D i V k dx + 0(R- 1 ) + 0(R- x ^ +3 ) 
Jb r * Jb r ijj>k 

= - 1) / h jDA jDM'l.r + o(l) 

as in Case 1. The rest of proof is the same as in Case 1. □ 
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